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The low-temperature crystal structure of BaCuSi2 06 has been investigated with high-resolution 
synchrotron x-ray and neutron powder diffraction techniques and has been found to be on aver- 
age (ignoring the incommensurate modulation) orthorhombic, with the most probable space group 
I bam. The Cu — Cu dimers in this material are forming two types of layers with distinctly differ- 
ent interatomic distances. Subtle changes also modify the partially frustrated interlayer Cu — Cu 
exchange paths. The present results corroborate the interpretation of low-temperature nuclear mag- 
netic resonance and inelastic neutron scattering data in terms of distinct dimer layers. 

PACS numbers: 61.05.C-, 61.05.fm, 61.66.Fn, 64.60.-i 



I. INTRODUCTION 



Intense studies of field-induced quantum phase tran- 
sitions (QPTs) in various magnetic insulators continue 
to enrich our knowledge of the possible quantum ground 
states of matter i~— Structurally dimer ized quantum spin 
systems play a leading role in these studies; the classes of 
field-induced QPTs known to date in these systems are 
well summarized in Refs. iand0 For magnetic interac- 
tions with weak or no frustration, the kinetic energy of 
the triplet quasiparticles in such systems is dominant and 
the ordered, or Bose-Einstein condensed (BEC), ground 
state is uniform at H > H C1 as e.g. in TlCuClsAi In 
SrCu2(BOs)2^^ triplet hopping is suppressed by geo- 
metrical frustration and the repulsion causes the con- 
densed triplets to form a superlattice with the appear- 
ance of magnetization plateaus. To better distinguish 
between the various scenarios and to determine accu- 
rately the relevant exchange couplings adequate knowl- 
edge of crystal structures and symmetries is of outmost 
relevance. 

The interest in BaCuSi2C>6 is motivated by its extraor- 
dinary phase diagram with field- induced BEC ^ 10 ^ 1 Be- 
ing a quantum paramagnet at zero magnetic field down to 
the lowest temperatures, the system displays a QPT into 
a magnetically ordered state at the critical value of mag- 
netic field of - 23.5 T 311 At - 610 K, the crystal struc- 
ture of BaCuSi206 undergoes a phase transition from the 
high-temperature (HT) phase with space group 14/mmm 
and lattice parameters aht ^ 7.11 A, Cht ~ 11.2 A to 
its room-temperature (RT) phase with a space group 
IAi/acd and lattice parameters a rt ~ y/2 * aht ~ 10 A, 
c r t ~ 2 * Cht ~ 22.5 A. 12 In the essentially layered RT 
tetragonal crystal structure of BaCuSi206 4 12lls the cop- 
per silicate Cu2(Si03)4 layers are separated by the in- 
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FIG. 1. (Color online) General view of the room tempera- 
ture tetragonal crystal structure of BaCuSi2 06- Copper sili- 
cate Cu2(SiOs)4 layers are interleaved with the layers of Ba 
atoms. The drawing is based on the crystal structure param- 
eters refined from neutron powder diffraction data collected 
at T = 130 K. 



termediate layers of Ba atoms (illustration in Fig. [I]). 
The Cu — Cu dimers with an interatomic distance of 
~ 2.75 A are well separated from each other in the struc- 
ture, since the inter-dimer Cu — Cu distances within the 
layer (~ 7 A) and the Cu — Cu distances between the 
dimers in adjacent layers (~ 5.75 A) are much larger. 
Thus, to a great extent, the ground state of the system 
is determined by the Cu — Cu interactions within the 
dimers . 14 i 15 

At zero magnetic field, the (S = 1/2) Cu 2+ ions are 
paired antiferromagnetically, while the degenerate ex- 
cited triplet states with 5 = 1 are gapped^ thus the 
compound possesses an essentially singlet ground state 
with S = at H = 0. Upon application of a mag- 
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netic field, the excited triplet states are Zeeman split 
with the energy scale of the splitting being proportional 
to magnetic field. At a critical magnetic field value of 
H c ~ 23.5 T, the energy of the lowest triplet state with 
S z = +1 becomes lower than that of a singlet state with 
5 = 0, and a magnetic field-induced BEC of the exci- 
tations occurs. The population of the bosons (measured 
as the magnetization of the sample) may be precisely 
tuned by the magnetic field. Another intriguing phe- 
nomenon observed in BaCuSi206 next to this QPT is 
dimensional crossover around 1 K from a 3D into a 2D 
regime with lowering temperature j 1 1 1 16 To explain this 
feature, perfect frustration between adjacent Cu2(Si03)4 
layers was assumed i 11 ! 16 ^ 7 Alternatively, the hypothe- 
sis has emerged that a structural modulation along the 
c axis^ i 18 i 19 combined with nonperfect frustration is the 
sole reason for the lowered dimensionality ! 20 1 2 1 As a uni- 
fying theory, a possibility has been proposed^ 2 - that the 
combined effect of the structural modulation along the 
c axis and the nearly perfect interlayer frustration is a 
reason for the observed dimensional crossover. 

In attempts to build a theory of this QPT, the com- 
pound was originally assumed to possess the same square 
motif of the Cu — Cu dimers at low temperatures as at 
room temperatur e) 10 ! 11 ! 23 resulting in three relevant cou- 
pling constants entering the spin Hamiltonian: the in- 
tradimer, inter dimer, and interlayer couplings. 

However, high resolution inelastic neutron scattering 
(INS) data were clearly at odds with this pictured Fur- 
thermore, a first-order structural phase transition at 
~ 100 K has been discovered by Stern et alJ& and re- 
ported by Samulon et al^ which necessarily states that 
the true spin Hamiltonian of the system is more com- 
plex than originally thought. The method of single crys- 
tal diffraction used in Ref. [24] confirmed the transition 
into an orthorhombic (or weakly monoclinic) structure. 
Moreover, it allowed for observation of satellite peaks 
indicating an incommensurate modulation in the low- 
temperature (LT) phase of BaCuSi206- Here we use 
in contrast a combination of high-resolution x-ray and 
neutron powder diffraction techniques in order to deter- 
mine the average crystal structure of the LT phase of 
BaCuSi206- 



II. EXPERIMENTAL 

A. Synthesis 

Polycrystalline samples of BaCuSi206 were prepared 
by solid-state reaction. Powders of BaC03, CuO, and 
Si02 were weighted to the prescribed ratios, mixed, and 
well ground. The mixture was calcined at 900° C in air for 
20 h. The resulting powders were pulverized, isostatically 
pressed into a rod shape (~ 5 mm diameter, ~ 50 mm 
length) and sintered again at 1010°C in air for 20 h. 



B. Powder diffraction experiments 

The powder diffraction experiments have been carried 
out with two high-resolution diffractometers: the Pow- 
der Diffraction station of the Materials Sciences Beam- 
line (MS-PD) at the Swiss Light Source, 25 and the high- 
resolution powder neutron diffractometer HRPT— at the 
spallation neutron source SINQ, both at Paul Scherrer 
Institute in Villigen. 

At the MS-PD station, the synchrotron x-ray diffrac- 
tion data were collected on a powder sample enclosed in 
a capillary of a 0.3 mm diameter, which was placed in a 
Janis flow- type cryostat. The Microstrip Mythen-II de- 
tector was used, which allowed for high counting rates 
while maintaining the high resolution which was essen- 
tially sample-conditioned. The typical counts of ^ 2* 10 5 
in the strongest peaks were achieved within ~ 1 minute. 
On cooling, the diffraction patterns were collected with 
1 K steps in the temperature range from 154 to 64 K, 
and with 2 K steps in the temperature range from 64 to 
12 K. After collecting the lowest temperature dataset at 
4 K, the data were collected on heating the sample in 
the temperature range from 5 to 123 K with typically 
1 ... 2 K steps and from 128 to 178 K with 5 K steps, in 
order to observe the transition both on cooling and on 
heating and to characterize the hysteresis thereof. 

At HRPT, a bigger amount (~ 1 g) of the very same 
powder sample of BaCuSi206 was enclosed into a vana- 
dium can of a 6 mm diameter, and the data were collected 
using a close-cycle refrigerator at 130 K (well above the 
transition on cooling) and at 13 K (far below the tran- 
sition) in the instrument setup with the highest resolu- 
tion. The ultimate resolution was needed because the 
orthorhombic splitting of the peaks of the parent tetrag- 
onal RT structure characterized by an orthorhombic- 
ity parameter defined as 2 * \a orth - b or th\/(a rth + b or th) 
amounts only to ~ 1.7 * 10 -3 . 

The data of the neutron diffraction experiment were 
used for the structure determination, and for the pre- 
cise refinement of the structural parameters at low tem- 
perature, while the synchrotron x-ray powder diffraction 
data have mainly been used for indexing, space group 
selection and refinement of the temperature dependence 
of the main structural parameters, as well as for iden- 
tifying all the impurity phases present. The following 
impurities have been determined in the powder sample: 
BaCu 2 Si 2 7 (~ 4.4 % wt.), BaCuSi 4 Oi (~ 1.3 % wt.), 
Ba 4 Si 6 0i 6 (~ 1.1 % wt.), BaSi0 3 (~ 0.3 % wt.), and 
CU2O (~ 0.1 % wt.). Approximately 93% wt. of the pow- 
der sample is the main BaCuSi206 phase. While treat- 
ing the neutron data with ~ 5000 counts in the strongest 
peak, the only relevant impurity which was definitely seen 
in the diffraction pattern and correspondingly input into 
the refinements was BaCu2Si2 07 (~ 4.4 % wt.); the oth- 
ers were undetectable, and thus were not considered. 

The RT crystal structure model (Ref. [T2I ) has been un- 
ambiguously confirmed by the Rietveld refinements based 
on both synchrotron x-ray and neutron data above the 
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FIG. 2. (Color online) Selected ranges of the synchrotron x- 
ray powder diffraction patterns of BaCuSi2 06 taken with 1 K 
steps on cooling from 120 to 64 K. Obvious are the jump in 
the unit cell volume at the transition at ~ 89.7 K (all the 
main peaks are shifted to higher 2G angles) and the splitting 
of the diffraction peaks with h / k. 

phase transition temperature. We do not focus on this 
issue here, but rather concentrate on the LT structure. 



III. RESULTS 

A. Low-temperature crystal structure of 
BaCuSi 2 6 

The LT structural phase transition in BaCuSi206 is 
clearly seen in both neutron and x-ray diffraction pat- 
terns, and is most evident in the temperature-dependent 
x-ray data. Illustrations of the transition are given in 
Figs. [2] and G3 While all the peaks get a little narrower 
when cooling below the transition, we do not observe 
any of them to broaden. Instead, some remain single 
peaks, and some are splitting. The peak (2,2,12) at 
~ 17.35° remains single, while the (6,0,4) and (6,2,0) 
at ~ 17.9° and ~ 18.15° become split into the or- 
thorhombic pairs (0, 6,4)/(6, 0,4) and (2, 6, 0)/(6, 2, 0), 
correspondingly. The peak (8,0,0) at ~ 23° splits into 
a (0,8,0)/(8,0,0) doublet, and the strong (6,4,8) at 
~ 23.1° transforms into a pair of (4, 6, 8)/(6, 4, 8) which 
are not resolved, since the amplitudes of their scattering 
vectors are too close to each other. The backward phase 
transition on heating manifests itself (see Fig. [3]) as a re- 
covery of the diffraction pattern to that of the RT crystal 
structure, yet it occurs at a higher temperature. 

The LT synchrotron x-ray diffraction pattern of 
BaCuSi206 showing clear splitting of the diffrac- 
tion peaks (hkl) with h ^ k of the parent room 
temperature tetragonal crystal structure (space group 
IAi/acd) has been indexed on the orthorhombic unit 
cell with lattice constants (at T = 4 K) {a, 6, c} = 
{9.951 A, 9.968 A, 22.239 A}. From the systematic ab- 
sences it has been found that the most probable space 
group is body-centered, and all the relevant possibilities 
have been checked for the structure determination. In 
fact we have performed the structure determination for 




FIG. 3. (Color online) Selected range of the synchrotron x- 
ray powder diffraction patterns of BaCuSi2 06 taken with 2 K 
steps on heating from 83 to 97 K and with 1 K steps on heating 
from 99 to 119 K (note the nonuniform scale on the temper- 
ature axis). One clearly sees the structural phase transition 
occurring (on heating) at ~ 107.2 K. 



all space groups and combinations of (a, 6, c), which did 
account for all diffraction peaks observed, although for 
many of these some of the calculated peaks were actually 
absent in the experimental pattern. Having done so, we 
may exclude any chance for having missed the true solu- 
tion. The structure determinations have been carried out 
on a neutron diffraction dataset taken at T = 13 K, either 
with the program FOX 27 or by direct symmetry reduc- 
tion in cases where it was possible and straightforward. 
All the models have been checked with extensive Rietveld 
refinements (done with the program FULLPROF 28 ) of their 
parameters from both the synchrotron x-ray and neutron 
powder diffraction patterns. 

While going from higher to lower symmetries in the 
orthorhombic syngony, already the space group Ibca 
(No. 73), which is a direct subgroup of the RT structure 
space group IAi/acd, gives the first sensible solution. Yet 
the next group Ibam (No. 72) supplies a definitely better 
model. Compared to the Ibca model, the model in Ibam 
provides a significant improvement of the agreement fac- 
tors: x 2 decreases from 4.33 to 2.66, and the Bragg R- 
factor - from 7.05 to 4.80. Both structure models are in 
fact slight distortions of the RT tetragonal structure with 
IAi/acd symmetry, but are having very essential differ- 
ence between each other from the point of view of the 
geometry of the Cu — Cu dimer layers and their stacking 
along the direction perpendicular to the Cu2(SiOs)4 lay- 
ers. The unique Cu atom in the RT structure model is 
being split into two distinct atoms in both models, but in 
the case of the Ibca model they are forming just one type 
of Cu— Cu dimer layer, while the correct Ibam model con- 
tains two distinctly different Cu — Cu dimer layers, and 
these will be discussed below. The Rietveld refinement 
plot obtained with the neutron data is shown in Fig. HI 
The refined parameters of the LT crystal structure are 
given in Table HJ and the most relevant interatomic dis- 
tances in Table [III The schematic representation of the 
crystal structure of the LT phase of BaCuSi206 is given 
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FIG. 4. (Color online) Rietveld refinement of the crystal 
structure parameters of the low-temperature (LT) phase of 
BaCuSi2 06 from the neutron powder diffraction data. The 
observed intensity, calculated profile, and difference curve 
are shown. The rows of ticks at the bottom correspond to 
the calculated diffraction peak positions of the phases (from 
top to bottom): BaCuSi2 06 - orthorhombic LT structure, 
BaCu 2 Si 2 07 impurity (~ 4.4 % wt.) 




FIG. 5. (Color online) Schematic representation of the crystal 
structure of the LT phase of BaCuSi206- In the right panel, 
only the Cu— Cu dimers are shown. The Cul — Cul distance of 
2.774(6) A is much longer than the Cu2-Cu2 with 2.701(6) A. 
For the inter-layer coupling paths, Cul — Cu2 (shown in thin 
lines in the right figure) , there are always two pairs of slightly 
inequivalent bond distances (see Table ITT)) . 



TABLE I. The refined crystal structure parameters of the 
LT phase of BaCuSi2 06- Results obtained from the refine- 
ment done on neutron powder diffraction data collected at 
13 K with the wavelength A = 1.8857 A. Space group Ibam 
(#72, Z=16). Isotropic temperature factors were refined with 
a constraint to equality for Cu atoms, for Si atoms, for the 
01... 04 atoms ("apical" ones of the Si04 tetrahedra, also 
coordinating Cu in the Cu04 squares), and for the 05 . . . 07 
atoms (directly bridging the Si atoms, i.e. located at exactly 
or roughly the same z values as Si). Atoms are located in the 
following crystallographic positions: Cul - in the Si (0, 0.5, z), 
Cu2 - in the 8h (0, 0, z) , Sil - in the 8g (0, y, 1 /4) , Si2 - in the 
8/ (x, 0,1/4), Si3, Si4, 06 and 07 - in the 8j(x,y,0) sites 
correspondingly; all the others - in the general 16k positions 
(x,y,z). All positional parameters were refined without any 
constraints. 



BaCuSi2 06 Lattice parameters at T = 13 K 
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TABLE II. The most relevant interatomic distances ( A) in 
the crystal structure of the LT phase of BaCuSi2 06 refined 
from the neutron powder diffraction data at T = 13 K. 
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in Fig. O 

Unlike in the single crystal experiments reported in 
Ref. [24], our study carried out by powder diffraction could 
not reveal the presence of any incommensurate peaks due 
to their relative weakness: According to Ref. \2% the cen- 
tral Bragg peaks did have to be scaled down by a factor 
of ~ 10 4 to display on the same scale with the satellite in- 
commensurate modulation peaks. Such dynamic ranges 
are inaccessible for powder diffraction, and in this sense, 



the presence of the incommensurate modulation did not 
hinder the determination of what we believe to be a most 
probable model for the average LT crystal structure of 
BaCuSi 2 6 . 

An interesting feature is observed in the geometry of 
the Cu2(Si03)4 layers which is very relevant for the inter- 
pretation of low-temperature magnetic properties: Mu- 
tual rotation of the CUO4 coordination polyhedra in the 
Cu-dimer layers, characteristic for the RT structure on 
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T=13 K 
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Cu1-Cu1 dimers 



FIG. 6. (Color online) Cu-Cu dimer layer structure, (a): In 
tetragonal phase (the figure is based on the atom positions 
refined from the neutron data at 130 K). Cu — Cu distance is 
2.731(2) A. (b): Two types of Cu — Cu dimer layers at low 
temperature as refined from the neutron data taken at 13 K. 
The Cu2-Cu2 distance is 2.701(6) A; the Cul-Cul distance 
is 2.774(6) A. 



BaCuSi2 06, is only preserved in the layers with shorter 
Cu — Cu dumbbell distances. This feature is illustrated in 
Fig. [6l At 130 K, in the tetragonal phase, all the Cu — Cu 
dimer layers are equivalent and the CUO4 coordination 
polyhedra look like squares slightly turned (~ 19°) with 
respect to each other in each Cu — Cu dumbbell around 
its axis [Fig. El (a): view along the c direction of the unit 
cell]. At low temperature, only the Cu2 — Cu2 dimer 
layer with shorter Cu — Cu distances, shown on the left 
of Fig. El (6), preserves this feature (the mutual rotation 
angles of the Cu204 squares with respect to each other 
are on the order of ~ 21 . . . 22°), and the Si40i2 rings in 
these layers are also (same as at 130 K) nearly squares. 
In contrast to this, in the Cul — Cul dimer layer with 
longer Cu — Cu distances, shown on the right of Fig. El (6), 
the CUO4 coordination polyhedra are not turned with re- 
spect to each other. The Si40i2 4-member rings in the 
Cul — Cul dimer layers show clear deviations from the 
ideal square shape, forming a checkerboard type of order- 
ing of their elongation directions along one of the a — b 
plane diagonals. The intralayer Cu — Cu distances are 
identical in both types of layers in our model, yet it is 
too hard to quantify the obvious differences in the exact 
intralayer exchange paths between the dimers, which are 
obviously nonequivalent for the layers with the shorter 
and the longer Cu — Cu dimers. 
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FIG. 7. (Color online) Unit cell volumes on cooling and on 
heating (a) and the relative phase weight proportions (b) for 
the tetragonal and orthorhombic modifications of BaCuSi2 06- 
Results of the refinements done with the synchrotron x-ray 
data. Symbol sizes (T) and (O) in panel (a) are proportional 
to the corresponding weight phase fractions. Error bars in 
(a) are also shown and are typically much smaller than the 
symbol sizes. 



B. First-order character of the transition 



In agreement with the results already reported 
in Ref. 0, the structural transition we observe in 
BaCuSi206 is of first order. The refined unit cell vol- 
umes, as well as the weight proportions of the tetragonal 
and orthorhombic phases on cooling and on heating, are 
shown in Fig. The observed discontinuity in the unit 
cell volume at the transition is ~ 1%. The 50 : 50 pro- 
portions of the phases are reached at ~ 89.7 K on cooling 
and at ~ 107.2 K on heating. The width of the hysteresis 
is thus ~ 17.5 K. 
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C. Possible alternative structure models 

Dealing with the neutron powder dataset with limited 
intensity does not in principle exclude the possibility of 
having missed the lowering of the symmetry just due to 
the fact that the very weak diffraction intensities may 
potentially be confused with systematic absences. Hav- 
ing carried out an exhaustive search in the orthorhombic 
syngony for alternative crystal structure models of the 
LT phase of BaCuSi206, and upon carrying out care- 
ful refinements of their parameters, we can as well list 
all other possible structure models providing comparable 
agreement factors. It turns out that all alternative possi- 
bilities which deserve being considered are limited to the 
subgroups of the space group Ibam. 

In Table IIII| we have summarized all hypothetically 
possible models for the structure of the LT phase of 
BaCuSi206 providing plausible refinements, along with 
the key values - the numbers of the independent unique 
Cu— Cu dimer layers, the lengths of Cu— Cu dimers, some 
statistical indicators for the complexities of the models, 
and the corresponding agreement factors of the refine- 
ments. We find it inappropriate to present here all the 
structure models in details. The group-subgroup trans- 
formations are straightforward, and any models will be 
supplied on request. The general trend is that all al- 
ternative possibilities provide a pattern of the Cu — Cu 
dimer layers either very similar to or even more com- 
plex than that in our model with the space group Ibam 
(Tables HI and HI]). For some of them we have not 2 but 
4 independent Cu — Cu dimer layers. Since the appar- 
ent improvement in agreement factors of the refinements 
for lower symmetries is only achieved when the number 
of independent atoms and refinable parameters is get- 
ting very high, and due to the obvious similarity to the 
highest symmetry model with the space group Ibam (Ta- 
bles U and |TTJ) , we strongly tend to believe that the true 
structure of the LT phase of BaCuSi206 may sufficiently 
well and adequately be described by our model with the 
space group Ibam. Some later studies with higher preci- 
sion may one day show slight deviations from this model 
to one of those listed in Table Hill 



IV. DISCUSSION AND SUMMARY 

We have investigated the low-temperature structural 
phase transition in BaCuSi206 by means of powder syn- 
chrotron x-ray and neutron diffraction and have deter- 
mined the most probable average crystal structure of its 
LT phase. The proposed crystal structure model of the 
LT phase of BaCuSi206 has space group Ibam, and con- 
tains two individual Cu atom positions (with only one in 
the room-temperature structure). These two Cu atoms 
are forming two types of Cu — Cu dimers thus altering the 
most essential magnetic Cu — Cu exchange path lengths 
in a stacking manner along the c axis. The interlayer 
Cu — Cu exchange distances are also becoming barely dif- 



TABLE III. List of all possible space groups for the crystal 
structure of the LT phase of BaCuSi2 06- SG and Set de- 
note the space group number and particular setting. N a and 
N p are the total numbers of independent atoms and of re- 
fined positional parameters (atomic coordinates). Ncu and 
Nd denote the numbers of independent Cu atoms and of the 
unique Cu — Cu dimer layers in each model correspondingly. 
Cu — Cu distances for different dimer layers are separated by 
slashes. The agreement factors y^^Rp correspond to the re- 
finements based on neutron powder diffraction data collected 
at T = 13 K. 



SG 


Set 


N ai N p 


N Cu ,N D 


Cu - Cu, A 


x 2 ,R P 


72 
60 
57 
56 


Ibam 

Pbcn^ 

PbcrrP_ 

Pccn^ 


14, 30 
21,59 
25,63 
22, 58 


2,2 
2,2 
2,2 
4,2 


2.701(6)/2.774(6) 
2.713(6)/2.768(6) 
2.705(6)/2.770(6) 
2.71(2)/2.78(3) 


2.66, 4.80 
2.32, 4.41 
2.55, 4.69 
2.42, 4.53 


55 


Pbam^ 


26,62 


4,4 


2.71(2)/2.70(2)/ 
2.86(2)/2.70(2) 


2.59, 4.75 


50 


Pban^ 


24, 56 


4,4 


2.71(2)/2.77(3)/ 
2.70(2)/2.77(3) 


2.36, 4.45 


49 


Pccrn ^ 


28,60 


4,2b, 


2.73(3) and 2.82(2)/ 
2.69(3) ana 2.72(3) 


2.46, 4.59 


46 
45 
33 


I2cm ^ 

Iba2^ 

P2 lC n^ 


25,62 
22, 57 
40, 119 


2,2 
4,2 
4,2 


2.705(6)/2.767(6) 
2.71(3)/2.78(3) 
2.71 (3)/2.77 (4) 


2.50, 4.61 
2.36, 4.46 
2.26, 4.32 


32 


Pba2^ 


44, 115 


8,4 


2.78(6)/2.76(8)/ 
2.67(8)/2.80(8) 


2.23, 4.31 


30 


P2an^ 


42,117 


4,4 


2.73(3)/2.77(3)/ 
2.68(3)/2.77(3) 


2.19, 4.23 


on 

z9 


PcaZi _ 


/in im 
40, 119 


A O 

4, 2 


z.74(o)/z.7o(o) 


O 1 O AIT 

z.lz, 4.17 


28 


P2cm ma 


50, 125 


4,2^ 


2.67(4) and 2.74(4)/ 
2.76(5) and 2.78(5) 


2.29, 4.34 


27 


Pcc2^ 


44, 115 


8,2^ 


2.76(9) and 2.78(9)/ 
2.76(6) and 2.69(6) 


2.14, 4.20 


26 


P2 1 am^ 


48, 127 


4,4 


2.69(3)/2.70(3)/ 
2.86(2)/2.70(3) 


2.36, 4.42 


23 


7222^ 


24, 56 


4,2^ 


2.70(3) and 2.71(3)/ 
2.75(3) 


2.43, 4.51 


18 


P22i2i^ 


42, 118 


4,2^ 


2.67(3) and 2.74(3)/ 
2.77(5) 


2.23, 4.28 


16 


P222 


48, 112 


8,4^ 


2.85(4) and 2.70(4)/ 
2.78(6) and 2.75(6)/ 
2.68(4) and 2.62(5)/ 
2.78(6) and 2.75(6) 


2.22, 4.25 



a Subgroup of space group 72 (Ibam). 

b Different Cu — Cu distances within one or more dimer layers. 
c Subgroup of a few space groups which are subgroups of space 
group 72 (Ibam). 



ferent at low temperature, forming in total two distinct 
but close interlayer Cu — Cu distances. In agreement 
with the previous study (Ref. H3), we confirm the phase 
transition to be of first order, with a unit cell volume 
discontinuity of ~ 1% at the transition temperature, and 
a temperature hysteresis of ~ 17.5 K (the 50 : 50 phase 
weight proportions being achieved at ~ 89.7 K on cooling 
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and at ~ 107.2 K on heating). In the LT crystal struc- 
ture, the layers with shorter Cu — Cu dimer distances 
(~ 2.701 A at 13 K) preserve the mutual rotation of the 
CUO4 coordination squares with respect to each other 
in the Cu — Cu dumbbells and the squarelike shapes of 
the (SiOs)4 molecules around them. This feature which 
is characteristic for the room-temperature structure is 
lost in the layers with longer Cu — Cu bonds (~ 2.774 A 
at 13 K); the coordination CuC^ squares in them unfold 
the mutual rotation and are aligned identically, and the 
(Si03)4 molecules get strongly elongated along one of the 
a — b diagonals, forming a checkerboard pattern of their 
elongation directions. Thus, in the LT crystal structure 
of BaCuSi206, the two distinct adjacent Cu2(Si03)4 lay- 
ers with different Cu — Cu dimer distances and even dif- 
ferent internal geometries are alternating in a sequential 
manner along the c axis. This finding corroborates the 
interpretation of the low-temperature INS data^ in terms 
of existence of two strongly inequivalent dimer layers in 
BaCuSi206 with possible further complex modifications 
to one of them. The lower energy excitations reported 
in£ would correspond to the Cu2 — Cu2 dimer layers with 
shorter Cu — Cu distances, while the higher energy mul- 
timode excitations correspond to the Cul — Cul dimer 
layers with longer Cu — Cu distances. The incommensu- 
rate modulation of the crystal structure observed before 
in different studies did not hinder the structure determi- 
nation, since as shown in Ref. |24j the intensities of the 
satellite reflections are ~ 10 4 times weaker than that of 
the fundamental reflections, and thus they did not dis- 
turb our powder diffraction based average structure de- 
termination. 

Here we want to conclude with a remark about the pos- 
sibility of Dzyaloshinskii-Moriya (DM) interactions for 
BaCuSi 2 6 . All three, HT, RT, and LT, structures of 
BaCuSi2 06 are centrosymmetric. However, at RT there 



is no inversion center in the middle of the Cu— Cu dimers, 
as correctly pointed out by Ref. [24]. The same holds for 
LT Cu2 — Cu2 dimer layers with shorter Cu — Cu dis- 
tances. In contrast to this, the LT Cul — Cul dimer 
layer with longer Cu — Cu distances has similarly to the 
HT IA/mmm structure the inversion symmetry in the 
middle of the Cu dimer and the intradimer DM coupling 
would be forbidden. In the real structure of BaCuSi206, 
the incommensurate modulation may well remove this 
local symmetry and make DM interaction allowed again. 

We note that the results presented here provide a struc- 
tural model for the compound which justifies the inter- 
pretation of previous INS and NMR data by a number 
of groups*^ in terms of two inequivalent dimer layers in 
BaCuSi206 with further complex modifications to one of 
the layers that may partially relieve frustration. Further- 
more, the detailed structure information provided here 
will be important input for the ongoing ah initio calcu- 
lations of the exchange interactions in this fascinating 
material. 

ACKNOWLEDGMENTS 

This work is based on the experiments carried out at 
the Swiss Light Source synchrotron and at the Swiss 
Spallation Neutron Source SINQ, Paul Scherrer Insti- 
tute, Villigen, Switzerland. We acknowledge discussions 
with C. Berthier, M. Horvatic, S. Kramer, E. Joon, 
J. Mesot, C. Batista, F. Mila, A. Tsirlin, F. Due, and 
K. Sparta. Our special thanks go to C. Riiegg for numer- 
ous inspiring discussions and his careful reading of the 
manuscript. Work in Tallinn was supported by the Esto- 
nian Ministry of Education and Research under Grants 
No. SF0690029s09 and No. SF0690034s09, and Esto- 
nian Science Foundation under Grants No. ETF8198 and 
No. ETF8440. 



* Denis. Cheptiakov@psi.ch 
t lRaivo.Stern@kbfi.eel 

1 I. Affleck, Physical Review B 43, 3215 (1991). 

2 T. Giamarchi and A. M. Tsvelik, Physical Review B 59, 
11398 (1999). 

3 S. Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge, 1999) p. 353. 

4 T. M. Rice, Science 298, 760 (2002). 

5 C. Riiegg, D. F. McMorrow, B. Normand, H. M. Ronnow, 
S. E. Sebastian, I. R. Fisher, C. D. Batista, S. N. Gvasaliya, 
C. Niedermayer, and J. Stahn, Phys. Rev. Lett. 98, 017202 
(2007). 

6 T. Giamarchi, C. Riiegg, and O. Tchernyshyov, 
|Nat. Phys. 4, 198 (2008)] 

7 C. Riiegg, N. Cavadini, A. Furrer, H. U. Gudel, K. Kramer, 
H. Mutka, A. K. Habicht, P. Vorderwisch, and A. Wildes, 
Nature (London) 423, 62 (2003). 

8 H. Kageyama, K. Yoshimura, R. Stern, N. Mushnikov, 
K. Onizuka, M. Kato, K. Kosuge, C. P. Slichter, T. Goto, 
and Y. Ueda, Phys. Rev. Lett. 82, 3168 (1999). 



9 K. Kodama, M. Takigawa, M. Horvatic, C. Berthier, 
H. Kageyama, Y. Ueda, S. Miyahara, F. Becca, and 
F. Mila, Science 298, 395 (2002). 

10 M. Jaime, V. F. Correa, N. Harrison, C. D. Batista, 
N. Kawashima, Y. Kazuma, G. A. Jorge, R. Stern, I. Hein- 
maa, S. A. Zvyagin, Y. Sasago, and K. Uchinokura, Phys. 
Rev. Lett. 93, 087203 (2004). 

11 S. E. Sebastian, N. Harrison, C. D. Batista, L. Balicas, 
M. Jaime, P. A. Sharma, N. Kawashima, and I. R. Fisher, 
Nature (London) 441, 617 (2006). 

12 K. M. Sparta and G. Roth, Acta Crystallographica Section 
B-Structural Science 60, 491 (2004). 

13 L. W. Finger, R. M. Hazen, and R. J. Hemley, American 
Mineralogist 74, 952 (1989). 

14 Y. Sasago, K. Uchinokura, A. Zheludev, and G. Shirane, 
|Phys. Rev. B 55, 8357 (1997) | 

15 S. A. Zvyagin, J. Wosnitza, J. Krzystek, R. Stern, 
M. Jaime, Y. Sasago, and K. Uchinokura, Phys. Rev. B 
73, 094446 (2006). 

16 C. D. Batista, J. Schmalian, N. Kawashima, P. Sengupta, 



8 



S. E. Sebastian, N. Harrison, M. Jaime, and I. R. Fisher, 
Phys. Rev. Lett. 98, 257201 (2007). 

J. Schmalian and C. D. Batista, Phys. Rev. B 77, 094406 
(2008). 

M. Horvatic, C. Berthier, F. Tedoldi, A. Comment, 
M. Sofin, M. Jansen, and R. Stern, Progress of Theoretical 
Physics Supplement 159, 106 (2005). 

S. Kramer, R. Stern, M. Horvatic, C. Bert hier, T. Kimura, 

and I. R. Fisher, [Ph^s. Rev. B 76, 100406 (2007)] 

O. Rosch and M. Vojta, Phys. Rev. B 76, 180401 (2007). 

O. Rosch and M. Vojta, Phys. Rev. B 76, 224408 (2007). 

N. Laflorencie and F. Mila, 

|Phys. Rev. Lett. 102, 060602 (200971 

S. Sebastian, P. Sharma, M. Jaime, N. Harrison, V. Correa, 



L. Balicas, N. Kawashima, C. Batista, and I. Fisher, Phys. 
Rev. B 72, 100404 (2005). 

E. Samulon, Z. Islam, S. Sebastian, P. Brooks, M. Mc- 
Court, J. Ilavsky, and I. Fisher, Phys. Rev. B 73, 100407 

(2006). 

http:/ /www.psi.ch/sls/ms , 

P. Fischer, G. Frey, M. Koch, M. Konnecke, V. Pom- 
jakushin, J. Schefer, R. Thut, N. Schlumpf, R. Biirge, 
U. Greuter, S. Bondt, and E. Berruyer, Physica B 276, 
146 (2000). 

V. Favre-Nicolin and R. Cerny, J. Appl. Crystallogr. 35, 
734 (2002). 

J. Rodnguez-Carvajal, Physica B 192, 55 (1993). 



